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ABSTRACT 


Analysis  on  a regional  basis  of  the  NORSAR  seismic  bulletin 
data  for  the  three  years  from  April  1972  to  March  1975  shows 
that  there  has  been  no  significant  temporal  change  in  the 
detectability  during  this  time  period.  The  detectability  as 
been  examined  using  the  least  squares  cumulative  method, 
least  squares  incremental  method  and  the  Lacoss-Kelly  incre- 
mental method.  The  results  are  fairly  consistent,  although, 
the  incremental  methods  as  expected  give  higher  thresholds, 
for  example,  for  the  Japan-Kuriles-Kamchatka  region  the 

three  methods  give,  for  the  90%  detectability,  values  at 
NORSAR  magnitudes  3.8,  3.9  and  4.0,  respectively, 
values  have  also  been  estimated  using  a weighted  least 
squares  method  and  the  Utsu  method,  and  the  median  of  t e 
values  from  the  five  techniques  are  for  all  regions  in  t e 
range  0.7  to  1.1,  with  some  regional  variations.  An  extensive 
< .parison  with  PDE  solutions  from  USGS  shows  that  NORSAR 
■ eierally  reports  somewhat  smaller  magnitudes.  The  event 
locations  as  published  by  NORSAR  since  the  beginning 
of  1973  are  generally  better  than  before.  For  the  bes 
regions  the  median  location  difference  between  NORSAR  and 
USGS  estimates  is  about  100  km,  and  very  few  regions  show 
any  significant  difference  between  the  average  distances  an 
azimuths  as  reported  by  the  two  institutions. 


1.  INTRODUCTION 

The  construction  of  the  large  seismic  arrays  LASA  and  NORSAR 
during  the  time  period  1965-1971  were  major  events  in  the 
development  of  observational  seismology.  The  data  from  these 
arrays  have  been  used  in  the  investigation  of  a wide  range 
of  seismological  problems,  even  though  the  main  purpose  was 
the  more  limited  one  of  finding  ways  and  means  to  detect, 
locate  and  identify  small  seismic  events.  During  the  first 
few  years  after  the  arrays  were  completed,  a large  number  of 
reports  and  papers  were  published  in  which  the  capabilities 
of  the  arrays  in  this  respect  were  investigated.  However, 
these  studies  were  mostly  based  on  data  collected  during  the 
more  initial  phases  of  operation,  so  that  their  results  were 
often  affected  by  1)  later  improvements  and  changes,  and  2) 
not  having  collected  enough  data.  Also,  different  methods 
have  been  used  in  the  evaluations,  all  these  being  the 
reasons  why  many  relatively  different  results  were  produced, 
especially  on  detection  capabilities. 

The  latest  and  most  complete  report  on  the  capability  of 
the  NORSAR  array  to  detect  and  locate  seismic  events  was 
published  by  Bungum  and  Husebye  (1974)  . In  that  study,  the 
detection  capability  was  investigated  using  data  for  the 
time  period  April  1972-March  1973,  while  for  the  location 
capability  data  covering  the  time  interval  February-November 
1972  were  used.  Data  for  two  more  years  are  now  available, 
during  which  time  the  operational  situation  has  been  fairly 
stable,  and  the  purpose  of  this  paper  is  to  re-evaluate  some 
of  the  results  from  Bungum  and  Husebye  (1974) . Also,  several 
methods  will  be  used  this  time  to  evaluate  detection  capability, 
and  special  attention  will  be  given  to  the  problem  of  esti- 
mating the  recurrence  parameter  b. 
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2.  THE  NORSAR  OPERATIONAL  SITUATION 

We  refer  to  Bungum  et  al.  (1971)  and  Bungum  ana  Husebye  (1974) 
for  a comprehensive  description  of  the  NORSAR  array,  its 
physical  configuration  and  instrumentation,  and  its  data 
analysis  procedures.  However,  since  changes  (albeit  minor) 
have  been  implemented  also  after  Bungum  and  Husebye  (1974) 
published  their  paper,  a short  review  of  the  operational 
situation  is  nevertheless  warranted,  and  for  the  purpose  of 
completeness  we  start  with  the  initial  NORSAR  system. 

The  first  seismic  data  at  NORSAR  was  recorded  in  September 
1969.  Named  Interim  NORSAR,  it  included  only  a few  sensors, 
for  very  limited  time  intervals,  and  was  terminated  in 
September  1970.  At  that  time,  the  first  recording  using  all 
132  short  period  and  66  long  period  seismometers  was  initi- 
ated, and  one  started  to  experiment  with  the  first  full  array 
beam  deployment  (called  AB  306) . The  data  recording  intervals 
were  still  very  limited.  A new  beam  deployment  (AB  310)  was 
implemented  in  April  1971,  at  which  time  relatively  regular 
data  recording  and  analysis  started.  The  whole  year  of  1971 
was  characterized  by  extensive  software  development,  at  the 
same  time  as  the  recording  'down  time'  was  reduced  to  about 
3%.  In  December  of  that  year,  another  beam  deployment  (AB 
401)  was  implemented,  and  in  January  1972,  one  could  intro- 
duce a new  and  much  improved  set  of  regional  time  delay  and 
location  corrections.  The  same  month  al  >so  saw  the  imple- 
mentation of  a more  high  frequent  on-line  band  pass  filter 
(1.2-3. 2 Hz)  which  significantly  improved  the  signal-to- 
noise  ratios.  Because  of  all  these  changes,  the  evaluation 
by  Bungum  and  Husebye  (1974)  started  with  the  data  from 
February  1972. 

A large  number  of  modifications  have  been  introduced  also 
since  then,  although  their  overall  effect  has  been  less 
noticeable.  In  July  1972,  the  so-called  incoherent  beam- 
forming was  implemented  (Ringdal  et  al.,  1975),  which  im- 
proved the  capability  of  the  array  to  detect  events  within 
30°  epicentral  distance.  In  November  1972,  new  regional  time 
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delay  and  location  corrections  were  introduced,  although  the 
changes  this  time  were  minor  and  mainly  effective  only  in 
areas  of  low  seismic  activity.  The  same  can  be  said  about 
the  new  beam  deployment  (AB  411)  implemented  in  February 
1973. 

Other  changes  which  took  place  during  the  evaluation  period 
of  this  report  were: 

regionalization  of  filter  setting  in 
Event  Processor  (EP)  (February  1973),’ 

regionalization  in  EP  of  the  slowness  estimation  procedure 
(beampacking  vs.  beamforming)  (March  1973) 

- more  efficient  screening  in  EP  against  local  explosions 
(which  represent  a serious  false  alarm  problem)  (February 

1974) ; 

- implementation  of  weighted  beamforming  (Christof fersson 
and  Husebye,  1974)  in  EP,  to  be  used  on  specifically 
selected  weak  detections  (December  1974),  in  the  same 
package,  a number  of  statistical  tests  to  aid  in  the 
signal  identification  (Fyen  et  al.,  1975)  were  imple- 
mented (December  1974); 

- implementation  of  a floating  processing  threshold  in  EP 
(Steinert  et  al.,  1975)  for  the  purpose  of  keeping  a 
constant  false  alarm  rate,  thus  countering  diurnal  varia- 
tions in  the  noise  varianse  (December  1974) . 
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3.  THE  NORSAR  DETECTIC  ABILITY 

The  time  interval  used  in  the  detectability  calculations  by 
Bungum  and  Husebye  (1974)  was  from  /ipril  1972  to  March  1973, 
during  which  time  about  5500  events  were  reported  bv  NORSAR. 

The  data  were  grouped  into  15  regions,  a regionalization 
which  has  been  adopted  also  for  this  paper  (Table  1) . The 
first  data  column  in  Table  1 shows  the  number  of  events  analyzed 
by  Bungum  and  Husebye  (1974)  in  each  of  the  regions J it  is 
seen  that  for  two  of  the  regions  the  number  is  less  than 
100,  which  naturally  makes  a detectability  estimate  quite 
uncertain.  The  number  of  events  for  the  next  two  years 
(Table  1)  shows  a considerable  variation  for  some  regions, 
the  reason  for  that  is  most  likely  not  variation  in  detect- 
ability but  rather  variation  in  seismicity,  which  emphasizes 
the  importance  of  using  sufficiently  long  time  intervals  in 
this  type  of  analysis. 

Cumulative  Event  Detection  Capability  of  the  Array 

Using  the  same  method  as  Bungum  and  Husebye  (1974),  we  computed 
the  50%  and  90%  cumulative  detectability  estimates  for  each 
of  the  three  consecutive  years  covered  by  this  report,  as  well 
as  the  results  for  the  three  years  combined,  and  the  results 
are  shown  in  Table  2.  The  method  is  briefly  as  follows: 

A straight  line  xs  fitted  (see  Fig.  1),  using  the  method  of 
least  squares,  through  the  part  of  the  (cumulative  or  incre- 
mental) log  frequency-magnitude  distribution  which  is  between 
two  magnitude  limits.  These  limits  are  defined  by  the  analyst 
so  as  to  cover  a reasonably  straight  portion  of  the  distribution. 
The  90%  and  the  50%  detectability  limits  are  then  defined  as 
the  magnitudes  at  which  the  actual  number  of  events  falls 
10%  and  50%,  respectively,  below  the  level  predicted  by  the 
extrapolation  of  the  calculated  frequency-magnitude  distri- 
bution towards  lower  magnitudes . 

Although  there  is  a better  performance  for  some  of  the 
regions  in  the  last  period  (1974/75),  the  trend  is  towards 
slightly  higher  thresholds.  This  is  probably  not  caused  by  a 


real  deteriorating  performance,  but  rather  indicates  that 
the  false  alarm  problem  has  been  better  handled  in  the  later 
years,  thus  giving  fewer  events  reported  at  low  magnitudes. 

The  main  impression  is,  however,  that  there  is  a remarkable 
stability  in  the  performance  over  the  three  years,  thus 
giving  a consistent  data  base  for  threshold  determination 
combining  all  three  years.  As  an  example  of  the  fact  that  the 
changing  seismicity  level  does  not  significantly  affect  the 
estimates  of  detection  performance,  we  can  look  at  the  results 
for  region  9 (Ryukuo-Philippines) . In  Table  1 it  is  seen  that 
the  number  of  events  reported  for  this  region  is  961  in 
1972/1973  and  466  in  1974/1975,  while  Table  2 shows  that  the 
detectability  estimates  for  these  two  periods  are  essentially 
unchanged. 

The  obtained  results  show,  as  also  demonstrated  by  Bungum  and 
Husebye  (1974),  that  the  best  results  are  achieved  for  regions 
5-8  which  cover  the  belt  through  the  Mediterranean,  Iran, 

Western  Russia,  Central  Asia  to  Southeast  Asia.  All  of 
these  regions  are  located  on  the  same  (Eurasian)  tectonic 
plate  as  NORSAR.  Within  90  degrees  distance  range  from 
NORSAR  then  follow  the  Aleutian-Alaska  region,  Japan-Kamchatka 
Western-Northern  America  and  the  Mid-Atlantic  Ridge,  while 
the  poorest  performance  is  obtained  for  the  Central  America 
and  the  Ryukuo-Philippines  regions.  As  pointed  out  by  Bungum 
and  Husebye  (1974),  the  regions  where  NORSAR  has  the  best 
detectability  performance  are  the  regions  where  the  average 
signal  frequency  is  high,  while  generally  the  regions  where 
NORSAR  has  a poorer  performance  are  the  regions  from  where 
lower  frequencies  are  observed. 

Incremental  Detection  Capability 

It  is  sometimes  desirable  to  specify  the  incremental  detection 
capability  of  a station,  i.e.,  the  probability  of  detecting 
an  event  of  a given  magnitude.  It  is  clear  that  the  incre- 
mental detection  thresholds  should  always  be  higher  than  the 
cumulative  ones.  Lacoss  (1971)  computed  the  expected  differences 


Table  1 

Number  of  events  reported  by  NORSAR  for  the  intervals 
April  1972-March  1973,  April  1973-March  1974,  April  1974-March  1975 

and  April  1972-March  1975. 
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Cumulative  50  and  90  per  cent  detection  thresholds  in  terms  of  NORSAR  m^  units  for  the  intervals  April  1972- 
March  1973,  April  1973-March  1974,  April  1974-March  1975  and  April  1972-March  1975.  A dash  indicates  that 
there  were  too  few  data  or  a too  difficult  distribution  for  detectability  estimates. 
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CUMULATIVE 


fl=  7.37*0.03 

b=  -o.ee  *0.01 


INCREMENTAL 

fl=  6.31  *0.10 
B=-0. 80  *0.02 

50?  = 3. 84 
90?  = 3.92 


Frequency-magnitude  distribution  for  region  14 
for  April  1972  - March  1975.  Cumulative  and 
incremental  detectability  levels  are  indicated. 
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between  the  two  cases  based  on  a Gaussian  detectability  curve. 

He  found  values  typically  in  the  range  0.1-0. 4 m^  units,  and 
with  a somewhat  larger  difference  for  50  per  cent  thresho.*.s 
than  at  the  90  per  cent  level. 

Several  techniques  have  been  devised  to  estimate  the  incre- 
mental detectability  of  a seismic  station.  For  example,  the 
method  described  previously  for  the  cumulative  case  can  be 
applied  to  the  incremental  case  as  well,  and  the  results  are 
presented  in  Table  3.  It  should  be  noted,  however,  that  the 
estimates  by  this  method  are  quite  unstable,  mostly  because 
of  the  difficulties  inherent  in  estimating  the  incremental 
seismicity  slope  by  least  squares.  A different  approach  has 
been  used  by  Lacoss  and  Kelly  (1969),  who  applied  a maximum- 
likelihood  technique  to  simultaneously  estimate  seismicity 
and  detection  thresholds.  Briefly,  their  method  assumes  that 
the  number  of  earthquakes  N exceeding  a given  magnitude  m 
is  specified  by 

log^Q  Nc  = a - b • m (D 

and  that  the  detection  probability  of  a station,  given  the 
event  magnitude  m at  the  station  itself,  is  of  the  form: 

P (Detect/m)  = ^ -)  (2) 

where  $ is  the  standard  cumulative  Gaussian  distribution  func- 
tion. The  four  unknown  parameters  a,  b,  y and  o may  then  be 
estimated  by  maximizing  the  likelihood  function: 

K a-bm.  nu-y 

L(a,b,y,o)  = exp(-N)  • II  *10  • <*>(-— — ) (3) 

i=l 

Here,  mi,...,mK  represent  the  measured  magnitudes  of  all 
detected  events,  and  N is  the  total  expected  number  of  detected 
events : 


log1Q  N = a - by  + b2 •o2/(2*log10  e) 


(4) 
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Table  4 presents  the  results  from  incremental  threshold  estima- 
tion using  the  technique  of  Lacoss-Kelly . It  is  instructive 
to  compare  these  results  to  those  obtained  previously.  We  see 
that,  for  the  regions  with  the  largest  number  of  events,  the 
two  incremental  techniques  agree  fairly  well,  and  generally 
give  higher  results  than  the  cumulative  method.  For  example, 
looking  at  region  10  ( Japan-Kuriles-Kamchatka)  for  1972-75, 
we  obtain,  for  50  and  90  per  cent  thresholds,  respectively: 

3.4  and  3.8  (cumulative),  3.7  and  3.9  (incremental,  least 
squares)  and  3.7  and  4.0  (Lacoss-Kelly). 

It  should  be  noted  that  the  estimation  results  quoted  in  this 
section  always  refer  to  detection  thresholds  in  terms  of 
NORSAR  magnitudes.  As  shown  by  Ringdal  (1975),  the  curve  (2) 
can  in  principle  easily  be  converted  to  a 'true'  detection 
curve  (i.e.,  the  probability  of  detection  as  a function  of  a 
hypothetical  true  magnitude),  if  one  assumes  that  NORSAR  mb 
is  normally  distributed  (m+b^,  a^)  for  any  given  event  of 
true  magnitude  m.  The  true  detection  curve  then  becomes: 


m-y+bN 

P (detect/true  m)  = $( ) 


(5) 


Assuming  bN=0,  oN=0.25,  it  follows  that  the  incremental  50 
per  cent  thresholds  remain  unchanged,  while  the  90  per  cent 
values  determined  by  the  Lacoss-Kelly  method  must  be  incre- 
mented by  about  0.15  m^  units. 


Same  as  Table  2 for  the  incremental  detection  thresholds 


Table  4 

Same  as  Table  2 for  the  incremental  detection  thresholds 
using  the  Lacoss  and  Kelly  (19  69)  method  of  estimation. 
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4.  CALCULATION  OF  THE  RECURRENCE  PARAMETER  b 

As  a by-product  of  the  detectability  computations  in  the 
preceding  section/  we  have  obtained  estimates  of  the  seis- 
micity recurrence  parameters  a and  b in  the  formula  (1) • In 
fact,  large  arrays  are  particularly  well  suited  to  investigate 
regional  seismicity,  since  they  combine  a good  location  capa- 
bility with  magnitude  estimates  that  are  consistent  over 
the  whole  magnitude  range.  In  contrast,  seismic  networks 
tend  to  overestimate  magnitudes  of  small  earthquakes 
(Herrin  and  Tucker,  1972),  thereby  i ‘.storting  the  frequency- 
magnitude  relationship.  Furthermore,  the  variance  of  array 
magnitude  estimates  relative  to  the  "true  magnitude"  does  not 
bias  the  estimation  of  b,  and  can  be  compensated  for  when 
estimating  the  a parameter  (Ringdal,  1975)  . 


In  this  section  we  summarize  the  estimation  results  of  the 
parameter  b of  (1).  In  addition  to  the  three  methods  dis- 
cussed in  the  preceding  section  for  detectability  estimates, 
we  have  applied  Utsu's  (1965)  maximum  likelihood  estimate 


m-in 

o 


Here,  perfect  detectability  is  assumed  down  to  a magnitude  mQ, 
and  in  is  the  average  of  all  observed  magnitudes  exceeding  mQ. 

We  also  attempt  to  modify  the  least  squares  procedure  for 
incremental  estimation  by  weighing  each  data  point  (m^,  log  N^) 
by  the  weight  factor  (i.e.,  the  number  of  observations  at 

the  magnitude  itk)  . 


The  results  from  the  five  methods  of  ebtimation  are  shown  in 
Table  5.  The  entire  data  set  1972-75  has  been  used,  with  the 
same  regionalization  as  before.  We  note  that  the  results  show 
occasional  large  variation  between  the  different  techniques, 
while  some  common  trends  are  also  noted.  As  an  example,  the  two 
incremental  least  squares  methods  (2  and  3)  generally  produce 
the  lowest  b- values.  This  was  to  be  expected,  since  the 
slope  b'  of  the  'incremental'  relationship 


1 


-v- 


is  not  necessarily  equal  to  b in  (1).  It  is  clear  that  b=b' 
if  (1)  is  valid  without  any  restriction  on  m,  while  if  one 
assumes  a certain  highest  possible  magnitude  value,  then  the 
cumulative  slope  will  tend  to  decrease  more  steeply  than  the 
incremental  slope. 

There  is  clearly  a high  amount  of  uncertainty  involved  in 
interpreting  the  data  of  Table  5.  However,  we  note  that  the 
median  value  for  all  regions  lies  between  0.7  and  1.1. 
Examining  more  closely  the  regions  within  the  teleseismic 
zone,  we  see  that  regions  1,  2 and  10  (covering  the  seismicity 
belt  from  Japan  to  Western  North  America)  all  have  b-values 
at  or  below  0.80,  while  the  Eurasian  regions  5 and  6 exceed 
1.0.  Thus,  there  are  differences  of  some  significance  in 
these  data,  in  spite  of  the  large  amount  of  smoothing  intro- 
duced by  our  large-scale  regionalization.  It  would  be  of 
considerable  interest  to  use  NORSAR  data  to  obtain  b-values 
at  a more  detailed  regional  level.  Also,  it  would  be  inter- 
esting to  discuss  more  closely  the  relative  merits  of  tha 
various  estimation  techniques.  Both  of  these  topics  are, 
however,  outside  the  scope  of  the  present  paper. 
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5.  COMPARISON  BETWEEN  NORSAR  AND  USGS  BODY  WAVE  MAGNITUDES 

Comparing  magnitudes  reported  by  an  array  like  NORSAR  with 
those  reported  by  a network  like  USGS  is  a hazardous  task 
for  several  reasons.  If  NORSAR  did  not  have  any  amplitude 
saturation  problems,  the  best  procedure  would  be  to  use 
only  very  high  m^  events,  say  above  111^=6.  This  is  because 
for  such  events  a large  number  of  stations  and  not  only  the 
best  ones  would  report  to  USGS,  thus  one  would  not  expect  a 
bias  in  the  USGS  magnitude.  Secondly,  for  such  events  NORSAR 
would  have  essentially  100  per  cent  detectability,  and  no 
bias  would  be  introduced  by  the  fact  that  one  would  select 
only  those  events  detected  by  NORSAR,  and  thereby  only  those 
with  particularly  high  NORSAR  amplitudes.  NORSAR  does, 
however,  have  severe  clipping  problems,  and  since  it  is 
dxfficult  to  define  an  interval  where  we  expect  the  magnitudes 
to  be  unbiased,  we  have  used  all  events  jointly  reported  by 
USGS  and  NORSAR  in  the  period  January  1973-March  1975.  In 
Table  6 is  listed  the  average  magnitude  difference  together 
with  standard  deviation  and  number  of  events  used.  As  seen 
from  the  table  there  is  for  all  regions  a positive  difference, 
i.e.,  NORSAR  seems  to  be  reporting  smaller  magnitudes.  As 
pointed  out  by  Husebye  et  al.  (1974),  such  differences 
cannot  directly  be  explained  as  being  caused  by  beamforming 
loss  because  the  skew  amplitude  distribution  gives  a positive 
bias  to  the  array  beam,  which  roughly  compensates  the  loss 
caused  by  incoherent  signals.  In  Table  6 is  also  listed  the 
results  of  a linear  least  squares  fit  through  the  NORSAR/USGS 
magnitude  distribution,  where  a maximum  likelihood  method 
has  been  used.  As  seen  from  Table  6 there  are  clear  regional 
variations  in  the  NORSAR/USGS  magnitude  relationship.  It 
should,  however,  be  stressed  that  we  have  no  control  over 
different  types  of  'bias'  effects  influencing  the  results  in 
Table  6 . 


Estimates  of  average  and  standard  deviation  of  the  difference  between  USGS  and  NORSAR  body  wave 
magnitudes  together  with  the  standard  errors  in  the  estimates.  The  two  last  columns  give  maximum 
likelihood  estimate  of  X and  K in  the  relation  m^ (USGS) =X+K*m^ (NORSAR) . 
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THE  NORSAR  LOCATION  CAPABILITY 

The  difference  in  epicenter  solutions  has  been  calculated 
for  the  events  reported  jointly  by  USGS  and  NORSAR  in  the 
period  January  1973  until  March  1975.  Fig.  2 presents  the 
results  for  region  14  (distance  limits  30-90  degrees  from 
NORSAR)  in  increments  of  50  km  location  difference.  Fig.  3 
shows  the  same  data  for  the  interval  February  to  November 
1972.  It  is  seen  that  the  distribution  in  Fig.  2 has  a much 
shorter  tail,  indicating  that  most  of  the  large  location 
errors  have  been  removed.  Since  the  distributions  are  skew, 
the  50  per  cent  (median)  and  the  90  per  cent  level  of  location 
differences  are  used  as  characterizing  parameters.  In  Table  7 
the  values  observed  for  the  period  February  to  November  1972 
(Bungum  and  Husebye,  1974)  are  listed  together  with  the 
values  observed  for  the  period  January  1973  until  March 
1975.  For  regions  5,  6,  and  7 there  is  some  tendency  towards 
higher  values  in  the  last  period,*  this  is  because  no  detailed 
analysis  is  performed  for  some  of  the  nearest  events,'  the 
event  is  instead  located  at  the  beam  location  where  it  was 
detected.  For  all  the  other  regions  there  is  a clear  tendency 
towards  a reduction  in  location  differences.  This  is  partly 
due  to  the  implementation  of  new  region  corrections  in 
November  1972,  and  partly  due  to  more  experienced  analysts. 

In  Table  8 the  location  difference  has  been  split  into  dis- 
tance and  azimuth  differences.  Also  here  are  the  averages  as 
well  as  the  standard  deviations  usually  substantially  reduced 
in  the  last  period  relative  to  the  first  one.  It  is,  for 
example,  seen  that  in  the  last  period  only  one  region  has  an 
azimuth  difference  which  in  average  is  more  than  one  degree, 
and  only  two  regions  have  distance  differences  with  an 
average  of  more  than  100  km.  It  is  noteworthy  that  some  of  this 
scatter  necessarily  is  introduced  by  the  errors  in  the  USGS 
estimates  themselves.  Based  on  these  bias  considerations  as  well 
as  on  independent  theoretical  estimates  (Shlien  and  Toksoz,  1973) 
it  can  be  concluded  that  the  location  precisions  reported  here 
for  NORSAR  are  now  at  the  performance  limit  for  an  array  that 
size. 
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30-90  degrees  from  NORSAR) 
1973  until  March  1975. 
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7 . CONCLUSIONS 

Based  on  analysis  of  NORSAR  seismic  bulletin  data  for  the 

time  period  April  1972  - March  1975  we  conclude  that: 

1)  The  detectability  of  the  NORSAR  array  has  not  changed 
significantly  with  time  during  the  three  years  analyzed. 

2)  For  teleseismic  (third  zone)  events,  the  detectability 
is  highest  for  the  regions  between  the  Middle  East  and 
Southeast  Asia,  including  Central  Asia,  and  poorest  for 
America  and  the  North  Atlantic  Ridge. 

3)  The  90%  incremental  detectability  thresholds  for  the  best 
regions  are  in  the  range  3. 7-3. 9 measured  on  the  NORSAR 
magnitude  scale,  while  for  the  whole  distance  range 
30°-90°  the  value  is  around  4.0. 

4)  The  recurrence  parameter  b (the  slope  of  the  frequency- 
magnitude  distribution)  is  found  to  be  for  all  the  regions 
in  the  range  0.7  to  1.1. 

5)  The  relation  between  NORSAR  and  USGS  reported  magnitudes 
is  dependent  both  on  seismic  region  and  magnitude,  and 
include  several  bias  effects.  Generally,  however,  NORSAR 
reports  smaller  magnitudes. 

6)  The  average  location  difference  between  NORSAR  and  USGS 
reported  magnitudes  is  best  for  the  regions  Japan- 
Kamchatka,  Aleutians-Alaska  and  Central  Asia,  where  the 
median  location  difference  is  100,  110  and  120  km, 
respectively.  For  the  whole  distance  range  30  -90  from 
NORSAR,  the  corresponding  number  is  130  km. 

7)  The  possible  bias  effects  in  the  NORSAR  location  estimates 
are  negligible,  and  no  further  significant  improvement 

of  the  location  precisions  are  now  considered  possible. 


23 


REFERENCES 

Bungum,  H.  s E.S.  Husebye  (1974).  Analysis  of  the  operations 
capabilities  for  detection  and  location  of  seismic 
events  at  NORSAR,  Bull.  Seism.  Soc . Am_._,  64,  637-656. 

Bungum,  H.,  E.S.  Husebye  & F.  Ringdal  (1971).  The  NORSAR  array 
and  preliminary  results  of  data  analysis,  Geophys^ 

R.  Astr.  Soc.,  25,  115-126. 

Christof fersson,  A.  4 E.S.  Husebye  (1974).  Least  squares  signal 
estimation  techniques  In  analysis  of  seismic  array  recor 
P-waves,  Geophys.  J.  R.  Astr.  Soc.,  38.'  525  552. 

Fyen,  J.,  E.S.  Husebye  & A.  Christof fersson  (1975).  Sto.  ical 
classification  of  weak  seismic  signals  and  noise  at  the 
NORSAR  array,  Geophys.  J.  R.  Astr.  Sqc^,  42,  529-546. 

Herrin,  E.  & W.  Tucker  (1972).  On  the  estimation  of  bodywave 

magnitude,  Technical  Report  to  AFOSR,  Dallas  Geophys.  Lab., 

SMU,  Dallas,  Texas. 

Husebye,  E.S.,  A.  Dahle  S K.A.  Berteussen  (1974).  Bias  analysis 
of  NORSAR-  and  ISC-reported  seismic  event  mfa  magnitudes. 


J.  Geophys.  Res.,  79,  2967-2978. 

Lacoss,  R.T . (1971).  Seismic  detection  and  discrimination  - 

some  statistical  considerations,  Proceedings,  Seminar 
on  Seismology  and  Seismic  Arrays,  NTNF/NORSAR , Kjeller, 

Norway . 

Lacoss,  R.T . , S E.J.  Kelly  (1969).  Estimation  of  seismicity 

and  network  detection  capability,  HIT  Lincoln  Laboratories 

Technical  Note  41. 

Ringdal,  F.  (1975).  On  the  estimation  of  seismic  detection 
thresholds,  Bull.  Seism.  Soc.  Am.,  65,  1631-1642. 

Ringdal,  F.,  E.S.  Husebye  * A.  Dahle  (1975).  P-wave  envelope 
representation  in  event  detection  using  array  data 
in:  Beauchamp,  K.G.,  Exploitation  of  seismograph  Networks, 
Noordhoff-Leiden,  The  Netherlands. 


24 


Shlien,  S.  & M.N.  Toksoz  (1973).  Automatic  event  detection 

capabilities  of  large  aperture  seismic  arrays,  Bull.  Seism. 
Soc . Am. , 63,  1276-1288. 

Steinert,  0.,  E.S.  Husebye  & H.  Gjtfystdal  (1975).  Noise  variance 
fluctuations  and  earthquake  detectability,  J . Geophys .,  41_, 
289-302. 

Utsu,  T.  (1965).  A method  for  determining  the  value  of  b in 
a formula  log  n — a-bM  showing  the  magnitude-frequency 
relation  for  earthquakes  (in  Japanese  with  English  abstract) , 
Geophys.  Bull.  Hokkaido  Univ . , 13 . 


